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Abstract. We investigate the entanglement produced by a multi-path interferometer that is composed 
of two symmetric multiports, with phase shifts applied to the output of the first multiport. Particular 
attention is paid to the case when we have a single photon entering the interferometer. For this situation 
we derive a simple condition that characterize the types of entanglement that one can generate. We then 
show how one can use the results from the single photon case to determine what kinds of multi-photon 
entangled states one can prepare using the interferometer. 
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1 Introduction 

Entanglement is perhaps the most important and enig- 
matic feature of quantum mechanics. It is at the heart 
of quantum computing [TJ[5J[3] an d has many other appli- 
cations within the field of quantum information [4,5,6 . 
The task of generating entangled states is thus of practi- 
cal as well as fundamental importance. A simple way of 
generating entanglement, within optics, is to use a bal- 
anced beam splitter [7]. A straightforward, but powerful, 
generalization of this would be to use a Mach-Zehnder in- 
terferometer [5] . Using this simple setup one can generate 
many different entangled two mode states. For instance, 
this setup allows one to generate all possible two mode, 
single photon entangled states, i.e all states of the form 
a|l>|0) + 6|0)|lE 

It would seem natural to try to generate multipartite 
entangled states using a similar linear optical setup. This 
can be achieved using two symmetric multiports 11,12, 
fT5lfl6 , 17.18.19. 20 and phase shifts. Interferometers such 
as this have been studied previously [T'lllTB'lllZfj] and have 
been realized experimentally [16]. The advantage of us- 
ing this setup over more general schemes, such as |14) . is 



that it can be realized straightforwardly using commer- 
cially available fiber couplers |16j . The previous work on 
multi-path interferometers have focused on showing that a 
specific state or class of states can be generated using the 
setup. There has, however, not been a systematic study of 
the states that these interferometers can produce. In this 
paper we will give a general treatment that determines all 
possible states that can be generated for single photon in- 
put. These results will then be used to show how one can 
classify the possible output states for two photons entering 
the apparatus. 

The results will allow us to answer a number of basic 
questions regarding multi-path interferometers. For exam- 
ple, the Mach-Zehnder interferometer allows one to gener- 
ate all possible single photon (two mode) entangled states, 
does this hold in the higher dimensional case? We will 
demonstrate that this is not the case. Similarly, if one has 
a multi-path interferometer, can the setup function as a 
lower dimensional multi-path interferometer? For exam- 
ple, by ignoring one of the paths, can we choose the phase 
shifts so that a three path interferometer functions as a 
two path interferometer? Finally, we will investigate what 
kinds of two photon entangled states one can generate. 



a Corresponding author emaihthomas. brougham@gmail.com 
1 It is sometimes claimed that such states are not truly en- 
tangled. The states can, however, be used to violate a suitable 
Bell inequality [9"lllU| and thus are entangled in the conven- 
tional sense. Furthermore, the question: 'what is the photon 
entangled with?', can be addressed by noting that it is the two 
modes of the electromagnetic field that are entangled, not the 
photon. 



2 Multiports and multi-path interferometers 

Beam splitters are ubiquitous within quantum optics and 
are an essential component of many quantum information 
protocols [3l l22~[l23ll24] . Because of this it is interesting to 
consider generalizations of balanced beam splitters to the 
case where we have d inputs and d outputs. Such devices 
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Fig. 1. A schematic diagram of a multi-path interferometer. 
The elements Ml and M2 are symmetric multiports, while PS1, 
PS2,..., PS d - 1 and PS d are phase shifts. 



are known as symmetric multiports. If a single photon 
enters one of the inputs of a symmetric multiport, then 
we have probability 1/d of finding it at a given output. 
Without loss of generality, we can describe the action of 
a symmetric multiport on a single photon input by using 
the discrete Fourier transform. This means that the action 
of the beam splitter on the input state \ip) = J2k &fc|lfc)) 
where \l k ) = |0) 8)fe |l)|0)® d - 1 - fe , will be -> |V>') = 



ELo^'I 1 ^' where 



( 2mjk\ 
exp ( — - — I b k . 



7T xp l 



Using symmetric multiports we can describe a gener- 
alization of the Mach-Zehnder interferometer, which has 
d inputs and d outputs. The interferometer will consist of 
two symmetric multiports, with phase shifts performed on 
the outputs of the first multiport. This can be represented 
mathematically by the following unitary operator 



U = FPF, 



(2) 



where F is the discrete Fourier Transform, which we de- 
fine here as F = ^ mn ex^[2mmn/d)/^/d\l m )0-n\-> an d 
P represents the phase shifts acting on each mode, hence 
P = Efc ^fc|lfc)(lfc|) with |Afc| = 1. Suppose that a sin- 
gle photon enters one of the inputs of the interferometer. 
Without loss of generality we shall assume that the pho- 
ton enters the first input, hence the input state is |lo). 
The output state will thus be U\1q) — \ip) = E m c m|lm), 
where 

1 \ - / 27rimn\ 

c m = -£^K exp I — - — I . (3) 

n ^ ' 

Equation ([3]) can be expressed in the compact form c = 
FX/ yd, where the vectors c and A each have compo- 
nents Ck and Afc respectively, and F has elements F m 
exp(2irimn/ d) / \fd. The results are not substantially chang< 
by having the photon enter the interferometer via a dif- 
ferent input port. In particular, the set of possible output 
states that one can generate is not dependant on which 
input port the photon enters. For example, if we generate 
the state \ip) with the photon entering via the m-th input 
and with phase shifts A, then we can also generate this 
state by having the photon enter the first input and using 
phase shifts that are a permutation of A. 



The question we now ask is: can a given state of the 
form \tp) = ^2 k Cfc|lfc) be generated by the unitary trans- 
form ([2]). when a single photon enters the first input? If 
the state can be made, then the phase shifts needed to 
produce the state are given by 



VdFh 



(4) 



The fact that |Afc| = 1, for k = 0, 1, ...,d — 1, implies that 
the state can be produced if and only if |(.Ftc)fc| = 1/y/d. 
It will be convenient to write this condition in a different 
form. By using equation (j4]) together with the fact that 



I A t 



1, we can obtain 



2Trik\n 



= for k = 0, 1,2, ...d-l, 



^2 C ™ C *n exp 

(5) 

where we have used the fact that ^ m |c m | 2 = 1. Multi- 
plying the left hand side of ([5]) by exp(2iripk / d) and then 
summing over k leads to the expression 



E { 2irik\n — m + p] 
c m c n exp I 

k,m,n(m^n) 



= o. 



(6) 



(1) The only cases when the summation exp(27rifc[r 



p]/d) 



will not be zero is when n — m+p = 0orn — m + p = d. 
This observation allows us to obtain the following new set 
of equations 



d-l 

^ , c m c m -p 
m=0 



0; forp= 1,2,..., d-l, 



(7) 



where the indices of each c term are modulo d, hence 
c x +d — c x . It is clear that a solution of §5§ will also be 
a solution of ([7]). With some thought one can easily es- 
tablish that the converse is true and thus equation ([7]) is 
equivalent to equation ([5])- 



3 Characterizing the entanglement generated 
from a d-port interferometer 

While equation (0 totally characterizes the possible out- 
put states, there are situations where it is desirable to 
have an alternative condition. This can be illustrated by 
the following example. Suppose one wanted to generate an 
entangled state with coefficients given by c, which is not 
a solution of fl7J. Let us further suppose that it is possible 
to generate a state with coefficients c', where \c'j\ = \cj\, 
j = 0, 1, d— 1. It is clear that we can transform this par- 
ticular state to the desired one by simply making appro- 
^^priate phase shifts to the output of the interferometer. One 
6 could thus say that these two entangled states are equiv- 
alent. The sense in which the states are equivalent can 
be made precise by considering the set of concurrences, 
{C m ,n}, between modes m and n |27pl . It has been shown 



In general, multipartite entanglement cannot be fully char- 
acterized using only the concurrences between each particle [251 
126] . The entanglement of a single photon state, however, can 
be fully characterized using the set of concurrences [21] . 
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that for the single photon case, the concurrence between 
two modes m and n reduces to |21j 

Cm, n = 2|c m ||c rl 

I- (8) 

The probability distribution for finding the photon at a 
particular output port, P m — |c m | 2 , is thus enough to 
characterise the entanglement between any two modes of 
the output state. This raises the question of whether we 
can design the phase shifts of the interferometer so that we 
can produce a state with a specific probability distribution 
P m . For d = 2, it is clear that we can create states with 
any desired output probability distribution, i.e. we can 
prepare states of the form a|10}+6|01), where |a| 2 can have 
any value between and 1. For multi-path interferometers 
with d outputs, we will find that there are limits to the 
possible output probability distributions. 

In order to find general conditions on the allowed out- 
put probability distributions, we will use a simple geomet- 
rical representation of equation ([7|l. One can obtain this 
geometrical picture of ([7J by plotting the terms c m c* m _ p 
on the complex plane. Each term c m c* m _ v will correspond 
to directed lines of length |c m ||c m _ p |. The condition that 
the sum of these terms must equal zero, will correspond 
to the directed lines forming a closed polygon in the com- 
plex plane. It is clear that for any polygon the length of 
one side must be less than the sum of the lengths of the 
other sides. This fact together with equation ([7]) lead to 
the following condition 

|c m ||c m _p| < ^ \ c k\\ck- P \, rn,p= 1,2,..., d-1. (9) 

One can compare a specific photon probability distribu- 
tion against the necessary condition ©. If the probability 
distribution violates the set of inequalities, then we cannot 
prepare a state where the probability of finding the photon 
in a particular mode is the same as the specified probabil- 
ity distribution. The concurrence between any two modes 
of the output state is given by ([5]) ; we can use this result 
to re-cast equation © in the following form 

C m , m - P < ^ Ck,k- P , m,p = 1,2,... ,d- 1, (10) 

where the indices of the concurrence are taken to be mod- 
ulo d, i.e. C m +d,n+d — Cm.n- The inequalities (fT0|) place 
bounds on how entangled two modes of the output state 
can be. 

While we have shown that both equations (0 and (fTUl) 
are necessary conditions, we have not addressed the ques- 
tion of whether they are also sufficient. It will be proved in 
section 0] that for d — 3, the conditions are also sufficient. 
For d > 3 the question is still open, however, numerical 
investigation suggest that equations © and (fTU|) are not 
sufficient for d = 4 and 5. 

Some consequences of equation (| 10[) will be examined. 
Suppose we want to construct a state where only two 
modes were entangled. This would correspond to the con- 
currence being zero between all but the two modes that 



are entangled. Let us denote the two entangled modes as a 
and b. We thus have that C a ,b — Cb, a > 0, while C m ^ n = 
for all other choices of m and n. When the left hand side 
of equation (fTU|) is not equal to C a ,b or Cb. a , then the in- 
equality (fTU)) is trivially satisfied. Consider the situation 
where the left hand side of (fit)]) equals C a ,b- Clearly the 
integer p must equal \a — b\. If the two modes a and b 
are to be entangled then C a .b > 0, which implies that 
J2k Ck,k+b-a > 0, where k ^ a. This condition can only 
be fulfilled when the summation Ck.k+b-a, k ^ a, con- 
tains the term Cb, a , hence Ck,k+b-a = Cb. a , which implies 
that [26— a] mod d = [a] mod d. If this last expression is 
true then we see that \(2b — a) — a\ = 2\b — a\ — nd, where 
n is a positive integer. The fact that a and b belong to the 
set of integers {0, 1, 2, d— 1} means that \b— a\ < d and 
thus 2\b — a | = d, which implies that d must be even. If 
d is not even then C a ,b < and the two modes will not 
be entangled. We thus have the following result. An in- 
terferometer of the type described in equation © cannot 
create a state where only two modes are entangled, when 
the interferometer has an odd number of input and output 
ports. 

For d even, it is possible to prepare states where only 
two modes are entangled. For example the state \tp) = 
i|c ||l ) + |c 6 | lfe), where b = a + d/2 and |c a | 2 -I- |c 6 | 2 = 1, 
can be produced. This can be easily verified using equation 
©■ 

4 Entanglement generation using a 3-port 
interferometer 

The set of inequalities © are necessary conditions upon 
us generating a state with a photon probability distribu- 
tion |c m | 2 . When our interferometer has three inputs and 
outputs, it can be shown that these conditions are also 
sufficient. We will now give a proof of this fact. The ap- 
proach that we shall use is constructive in that it leads to 
a method for calculating the phase shifts need to engineer 
the desired state. 

When our interferometer has 3 input and output ports, 
equation ([7]) is a set of two equations, each of which is the 
complex conjugate of the other. Equation ([7]) thus reduces 
to a single equation 

CQC\ + ClCp + c 2 c* = 0. (11) 

The same argument that leads to equation © will yield 
the following set of inequalities 

M|c 2 | < |co||ci| + |ci||e 2 |, 
I cx | j cq 1 < |co||c 2 | + |ci||c 2 |, 
|ci||c 2 | < I c 1 1 C2 I + |c 1 |ci I . (12) 

If the inequalities (fl"2")l are not satisfied, then we cannot 
prepare a state with coefficients {c TO }. Suppose that we 
have a state with coefficients {c m }, which satisfy the in- 
equalities (fT2|) . but not equation (fTTj) . In this case plotting 
the complex numbers cqc^, c±Cq and c 2 c*, on the com- 
plex plane will not give a closed shape. The fact that the 
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lengths of the three directed lines satisfy the triangle in- 
equalities, equation (|12l) . means that they could be rotated 
so that they form a triangle. This is equivalent to chang- 
ing the phases of the complex numbers c m . We will now 
show that when d = 3, we can always choose the phases 
so that a set {|c m |} that satisfies (|12p will be solutions of 
equation (ITT1) . 

Let <j> m be the phase of c m , i.e. c m = \c m \e l ^ m and let 
7i = <Pa — 4*2 , 72 = 4>i — <Po and 73 = 4>2 — 4>i ■ The condition 
that a set of coefficients {c m } satisfy equation can be 
represented geometrically as saying that they must form 
a triangle in the complex plane, as shown in figure [5] The 
values of the angles a, b and c (that are defined in figure 
2) can be found using the cosine rule 



cos(a) = 
cos(6) = 
cos(c) = 



M 2 


C2 | 2 


+ 


|ci| 2 M 2 - 


|ci| 2 | 


col 2 
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co||ci||c 2 | 2 






m 2 : 
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|ci| 2 | 
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|ci| 2 M 2 - 
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|c 2 | 2 



(13) 



Using some simple geometry we find that 

a = 71 - 73 + 7T, 
b = 72 - 7i - i" 5 

C = 73 - 72 + 7T. (14) 

For the sake of simplicity we shall take </>o = 0, this is 
equivalent to making a global phase change on the output 
state = *^2 m c m \lm) ■ It is now straightforward to see 
that we can take the other phases to be 

4>i = 3 (26 + a + tt) , 
02 - i(6- a + 27r). (15) 
This is equivalent to preparing the state 

\1>) = M|lo) + |c 1 |xe^ +2fc )/ 3 |l 1 ) + |c 2 | a; 2 e I ( b - a )/ 3 |l 2 ), 

(16) 

where x = exp(iri/3). 

The problem of determining whether it is possible to 
engineer a state with a particular probability for finding 
the photon in a given mode, has been solved for the case 
of d = 3. The inequalities (fT2"j) represent necessary and 
sufficient conditions for one to engineer a state with the 
given probability distribution. 

Equation (Ti"2"|) tells us whether we can engineer a state 
with a particular property; this still leaves open the prob- 
lem of how to actually prepare the state, i.e. what phase 
shifts are required. This problem can be resolve by us- 
ing equation (|T^)) together with equation (jl]) to obtain the 
following phase shifts 

A = |c | + \ Cl \xe^ + ^/ 3 + \c 2 \x 2 e^'\ 

Ai = |c | + | Cl |xV( Q+2b >/ 3 + |c 2 |(* 2 )VC 6 - a )/ 3 , 

A 2 = |c | - | Cl |e l ( a+2b) / 3 + | C2 |e l ( fc - a >/ 3 , (17) 




Fig. 2. A diagram showing the terms C0C2, ciCq and C2c\ rep- 
resented as lines in the complex plane, where the lines form a 
triangle. 

where a and b are determined from equation (II 3[) . It is 
clear that this choice of phase shifts is not unique and 
that different phase shifts could be used to generate the 
same output photon probability distribution. 



5 Multi-photon entanglement generation 

From the perspective of quantum information, multi-photon 
entangled states are of great interest. For this reason we 
will investigate the types of multi-photon entanglement 
states that one can generate using the setup described in 
figure 1. 

Let a m and b m be bosonic annihilation operators that 
acts on the m-th mode, i.e. a m |0) = |l m ), 6 m |0) = |l m ), 

[aj,aj] = dijl and [&,„,&}>] = <WL The operators {a m } 
will be associated with the input modes of the field, while 
{b m } are associated with the output modes of the field. 
The fact that symmetric multiports are linear optical de- 
vices means that output operators {b m } are related to the 
output operators {a m } by the linear equations 

b m = yjKrafcOfc) (18) 
k 

where V m k are the elements of a unitary matrix, i.e. J2k Vmk 
Vnk = fimn- In order for us to make a connection be- 
tween V and equation @, we will consider the situation 
where a single photon enters the interferometer. The input 
state will be a^JO), while the output state is J2k ^mfc^lO), 
where we have used equation (I18[) . When we have a single 
photon input, we can use equation ^ to express the out- 
put state as U\l m ) = J2k ^m^jjO), hence V mk = U km i.e. 
U = V T . Equation ^ shows that U m o = (J2k Fmk^k)/Vd = 
c rn . For n ^ 0, a straightforward calculation shows that 
U mn = Cm+ni where we take the index of c to be modulo 
d, i.e. c x+ d — c x . Combining this with the previous results 
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we find that 



V m 



The single photon results, such as equation ©, can be 
used to place restrictions on the types of multi-photon 
entanglement that we can generate. 

We will illustrate the previous theory by investigating 
the entanglement that one can generate from a multi-path 
interferometer with three inputs and outputs. The first sit- 
uation that we will examine is when we have two photons 
that both enter through the same input. Without loss of 
generality, we can assume that the photons enter through 
input 0, i.e. our input state is |200) = (al) 2 \Q) /\/2. Using 
equation (1191) together with some simple algebra we find 
the following general form for the output states 



\ip) = N (eg 1 200) + c 2 |020) + c||002) 
+2c ci|110) + 2c c 2 |101) + 2c lC2 |011)), 



where TV is a normalization constant. Suppose we want to 
determine whether or not we can generate a two photon 
entangled state that has a particular probability distribu- 
tion for finding the photons in the outputs. For example, 
the state = a |200) + c*i|020) + a 2 |002) + a 3 |110) + 
O!4|101) + CK5IOII), has the probability distribution 



We can also consider generating entangled state by 

(19) having photons enter through two different inputs. With- 
out loss of generality, we assume that the photons enter 
through inputs and 1, i.e. the input state is 1 1 10) = 
agaj|0). From equation (fTU)l we find that the general form 
for the output state is 

= iv[c ci|200) +clc 2 |020) +c2c0|002) 

+(c c 2 + cl 2 )|110) + (cic 2 + c 2 )|101) + (coci + c|)|011) 

(23) 

where N is a normalization constant. We can again use 
our previous results to determine whether a given en- 
tangled state can be generated from this setup. Suppose 
that we again wish to generate an entangled state that 

(20) nas the probability distribution (|2"Tj) . From equation (1231) 
it is clear that |ao| = iV|co||ci|, |oi | — -/V|ci||c 2 | and 
|a 2 | = AT I Co 1 1 C2 1 - The triangle inequalities (|12l) can again 
be applied to determine whether a given state is allowed. 
The problem of determining what kind of three photon en- 
tangled states one can generate using the multiport, can 
in principle be tackled in the same manner. 



P200 : 
-Pi 10 



l^ol' 
led 



P020 : 
P101 



M 



Poo: 
P01 



\a 2 \ 

l"5| 



(21) 



6 Conclusions 



If a state with this probability distribution is to be gen- 
erated by having both photons enter a single input, then 
from equation (l20|) it is clear that \ck\ = y/\ak\/N, for 
k = 0,1,2, which implies that 



\a 3 \ = 2>/|ao||ai|, |a 4 | = 2 N /|a ||a 2 | & |a 5 | = 2i/|ai||a 2 



If these three necessary conditions are not satisfied, then 
the entangled state cannot be generated by having two 
photons enter the same input of a three port symmet- 
ric interferometer. When the conditions (|22p are satisfied, 
then we must check that the values of |cfe|, defined by 
\ck\ — y/\a.k\/N, satisfy the set of inequalities (fT"2"|) . If this 
last condition is met, then we will be able to generate an 
entangled state with the given probability distribution. 
The phase shifts need to generate the state will be the 
same as the phase shifts that generate the single photon 
state J2k c fc|lfc)i which are given in equation (fT7|) . 

The utility of the previous results are best illustrated 
by two quick examples. The output photon probability 
P200 = P020 = P002 = 1/9 and Pno = P101 = P011 = 2/9, 
will violate the conditions (|2"2"]l and thus we cannot gener- 
ate a state that have this output probability distribution. 
By contrast, the output probability distribution P200 = 
P020 = P002 = 1/15 and Puo = >ioi = Poll = 4/15 
will satisfy both equation (f2"2"j) and the inequalities (fT2"j) . 
Using equations (|13l) and (|17[) together with some simple 
algebra, we find that the phase shifts needed to obtain an 
entangled two photon state with the desired probability 
distribution are X = e™l 2 and Ai = A 2 = e~ 47r / 6 . 



We have studied the states that one can generate using 
a multi-path interferometer that is composed of two sym- 
metric multiports with phase shifts applied to the outputs 
of the first multiport. For a fixed input one can obtain 
different output states by changing the phase shifts. We 
gave a thorough treatment of the case when a single pho- 
ton enters the apparatus before investigating the case of 
(22)two input photons. 

For interferometers with three modes we where able to 
fully determine all the single photon states that one could 
generate in terms of a simple inequality on the concurren- 
cies between each mode. Furthermore, we determined the 
phase shifts required to generate a three mode entangled 
state with a prescribed output photon probability distri- 
bution. This is significant as it means the results can be 
applied directly to experimental realization of three path 
interferometers such as [16] . 

The results for the single photon case, while interesting 
in their own right, also have relevance in determining what 
types of multi-photon entanglement that one can prepare. 
We were able to use the single photon results to determine 
what types of two photon states could be prepared. 

Our results allowed us to address an number of basic 
questions relating to multi-path interferometers. The first 
of these was whether by adjusting the phase shifts one can 
prepare a single photon entangled state = J2 m c m\^m) , 
such that |c m | 2 realizes all possible d outcome probability 
distributions. For d — 2 (i.e. the Mach-Zehnder interfer- 
ometer) this is indeed the case, however, we have shown 
that this does not hold for d > 2. 

Another interesting finding was that in general one 
cannot use a higher dimensional multiport to simulate a 
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lower dimensional one. The significance of this is that one 12. 
cannot use a single experimental setup to perform a va- 
riety of interference experiments. We established this fact 
by showing that a three port interferometer (or indeed any 
interferometer with d odd) could not be used to simulate i3 - 
a Mach-Zchndcr interferometer. Interestingly, a four port 
interferometer can be used to simulate a Mach-Zehnder 
interferometer, but it cannot simulate a three port inter- 
ferometer. It may be the case that a d port interferometer ' 
can never simulate a d—1 port interferometer, while being ^ 
able to simulate some other lower dimensional interferom- 
eters. This question is, however, still unresolved. ^ 

In all of the cases we have studied the input to the in- 
terferometer was a Fock state. There are situations, such j~ 
as continuous variable encoding, where we could have co- 
herent states as inputs. It has been shown [25] that the 
transformation matrix for the input and output ampli- ig 
tudcs of the coherent states is the related to the trans- 20. 
formation matrix between the annihilation operators of 
the input and output modes. One can thus directly apply 
the results found in sections [3] and |4] to determine what 21. 
(unentanglecQ) output states can be generated for a given 
coherent state input. 22. 
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